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Abstract 


In the first part of this work a rationalisation and a generalisation of the approach to the analysis of the diffusive resistances in laboratory-scale fuel 
cells and electrodes was recommended and some interpretation tools for the purpose were suggested and some simple reference cases discussed. 

This second part is mainly devoted to some applications to experimental data. In particular, some sets of experimental data from Ansaldo Fuel 
Cells laboratories are examined. The data regard molten carbonate fuel cells at a 10cm x 10 cm laboratory scale; a first set of data was obtained 
under constant current and changing utilisation factor conditions; the second under limit current conditions. 

The analysis indicated that the cathodic diffusive resistances are not important and the anodic ones, when correctly determined, do not depend 
on the anodic gas flow rate, but are likely to be localised in the liquid phase. 
© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction e Given local kinetics of this kind, for instance, a linear irre- 
versible one 
Many fuel cell systems at the laboratory scale are large Ca j 
enough to require a more detailed analysis. Often, while the J = jso ( ) (: — +) exp(67o) (1) 
uniform temperature and voltage assumptions can be considered Cao 
appropriate, the uniform concentration assumption may corre- the statements 
spond to a gross over simplification. The analysis of a set of . 
laboratory data obtained from systems of this kind requires inter- ko = ( Jeo ) exp(Bno) (2) 
pretation instruments of an intermediate level, which, on the one FCao 
hand, overcome the simple electrode concept and, on the other j 
hand, avoid the cumbersome and time-consuming structure of r=- F" koCas = ke(Ca — Cas) = KCA (3) 
detailed large-scale models. 
In the first part of this work [1] some principal concepts were 1 i= 
presented and these can be summarised as follows: = (= iE =) (4) 


make it possible to reduce the electrochemical kinetics to the 
classical formulation of a first-order heterogeneous reaction 
at the boundary. 

e When the local kinetics is averaged on the surface of a finite 
electrode, the effect of the concentration field can be taken 
into account by the coefficient 


1 C Ca dS 
* Corresponding author. = An = Js , A>l (5) 
E-mail address: bosio@diam.unige.it (B. Bosio). A Cao CaoS 


e The local electrochemical kinetics of a finite electrode can be 
expressed in terms of an external, constant reference condition 
that yields a local kinetic expression that is often a function 
of a unique variable, for instance, the concentration or the 
utilisation degree of a key component of the fluid phase. 


0378-7753/$ — see front matter © 2007 Elsevier B.V. All rights reserved. 
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Nomenclature 


A 
Co 


C 
Ci 


vas 


Pe 


MS MA KTS RNY YN Ra ATS 


corrective coefficient, from the averaging opera- 
tion 

asymptotic (low current) corrective coefficient; 
see Eq. (B5) in Part 1 

see Eq. (19) 

volumetric concentration of the component i 
[kmol m~?] 

diffusivity [m? s~!] 

molar flow rate [mol s~!] 

current density [A m2] 

mean electrode current density [A m2] 
apparent kinetic constant [ms~!] 

mass transfer coefficient [m s7 !] 

flow dynamic coefficient [m s7!] 

intrinsic kinetic constant [m s7 !] 

global kinetic constant, local value [m s7!] 
length [m] 

gas-liquid partition coefficient 

molecular weight [kg mol~'] 

stoichiometric coefficient for electrons, absolute 
value 

stoichiometric coefficients 

molar flux [mol m~? s71] 

absolute pressure [N m~?] 

Peclet number 

volumetric flow rate [m? s7 !] 

reaction rate, per unit surface [kmol m~? s7!] 
limit reaction rate [kmol m7? s7!] 

mean reaction rate [kmo] m7? s7!] 

gas constant [J kmol~! K~!] 

Reynolds number 

diffusive equivalent thickness [m] 

electrode surface [m7] 

Schmidt number 

Sherwood number 

absolute temperature [K] 

utilisation factor 

velocity [m s7!] 

voltage [W AW] 

advancement degree [kmol m3] 

see Eq. (19) 

mol fraction of the key reagent 

see Eq. (19) 


Greek letters 


ESZ NDR 


reaction order 

symmetry factor 

apparent reaction order 

adimensional potential loss 

dynamic viscosity [kg m7! s7!] 

mol number variation in respect to the key 
stoichiometric coefficient 

density [kg m7?] 

adimensional axial coordinate 


Subscripts and superscripts 


an anode 

A key reagent (hydrogen at the anode, oxygen at the 
cathode) 

c concentration 

ca cathode 

cp central point 

CP condensed phases 

eq equilibrium 

G gas phase 

i generic component 

L limit 

LP liquid phase 

m mean value 

M maximum value 

o inlet 

T reference 

s electrode surface 

T total 

a exchange 

/ 


refers to the other electrode 


so that, for linear kinetics, the apparent constant ką and the 
limit constant ky, 


Tm K 


= = Tl ke 
Cay A 


ka Sao 
Cao AL 


i ky 


(6) 


do not coincide with the global constant K and the trans- 
port coefficient kę, respectively, and their relationship with 
the intrinsic reaction constant ko 


1 1 1 


= 7 
Aka ko T ALkL 0) 


becomes more complex. 
In particular, the coefficients A and Ay, also depend on a flow 
dynamic constant, the fluid flow rate per unit electrode surface 


(8) 


that is the maximum apparent constant compatible with the 
reagent availability. 

The behaviour of the electrode can be directly interpreted 
in terms of Eq. (4) only if it is “simple”, that is subject to a 
uniform composition field. Otherwise, what is observed is the 
result of at least two factors: the local diffusive resistances set 
and the concentration distribution on the electrode. The first 
effect can be properly determined only after estimating the 
second with an adequate physical—-mathematical description. 
In many cases a description of sufficient complexity can be 
obtained by using the classical instruments of the chemi- 
cal reactor theory. Appendices D and E are devoted to some 
relevant results of this approach. 
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2. Analysis and correlation methods 


The characterisation of a fluid flow electrode or cell usually 
consists of the measurement of the current-voltage characteristic 
curves. Each of these curves is obtained by keeping the cell 
geometry and the inlet fluid composition and flow rate constant. 


2.1. Constant flow analysis 


In most cases, the flow dynamic coefficient kf is known and 
constant, while the apparent constant (ka, that is the current) and 
the intrinsic constant (Ko, that is the voltage), as well as the limit 
constant (kz, that is the limit current) need to be measured. In 
the case of a simple electrode the problem ends here, but it is 
often advisable to carry the investigation a little further. 

In fact, in the first place, Eq. (4) may be inadequate to describe 
complex electrodes, in terms of geometry, flow pattern or, par- 
ticularly, non-uniform composition fields. In the second place, 
Eq. (4) is certainly wrong when the reagent availability begins 
to manifest its limiting effects (low kf, that is high utilisation 
factors). In the third place, also for low utilisations, Eq. (4) and 
its asymptotic solution may not be fully adequate and require the 
introduction of a corrective coefficient co different from unity 
(see Part 1, Appendix B). 

When Eq. (4) is doubtful, the reference to a formulation witha 
wider range of validity, such as Eq. (7), is useful. The differences 
between Eqs. (4) and (7) are often appreciable, even in the low 
current range, where the local linearisation in Eqs. (5) and (6) 
of Part 1 can be transformed in terms of mean density current 


1 J Cot 
(PE > 


foc ¥ — = (10) 


by allowing for the asymptotic correction Co at low current. In 
such a way, the concentration polarisation ne is still considered 
proportional to the mean current density J and inversely pro- 
portional to the limit value Jy, of the mean current density (that 
is proportional to the ratio k,/k,, between the apparent and the 
limit constant), but the proportionality coefficient is determined 
as Co/ and only a careful examination of the characteristic curve 
can suggest whether it is correct or not to assume a unitary value 
for the asymptotic constant co and, in that case, how the correct 
value differs from unity for the constant. Obviously, when the 
analysis of the concentration polarisations is made in terms of 
ratios, all this discussion can be avoided, but it is important to 
have a correct estimation of the absolute values of polarisation. 
In particular, if the experimental value of the ratio J/ņoc is identi- 
fied with JL, this leads to an overestimation of the limit current; 
on the other hand, if the polarisation is estimated from an inde- 
pendent measurement of the limit current, such as noc = J/BJL, 
this leads to its underestimation. 

Finally, once a correct value for kg is determined, the local 
transport coefficient ke can be determined in its turn only by 
correcting kg by means of a coefficient AL, usually different 
from (greater than) unity. 


2.2. The effects of the flow conditions 


In all cases, the structure of the limit constant ky, or, better, 
the transfer coefficient kę and the effects of the flow conditions 
on them is to be carefully examined, even if Eq. (4) stands. This 
can be done in an empirical way, by correlating experimental 
data on the limit current at different flow rates or different inlet 
compositions in terms of adimensional coefficients, for instance, 
by expressing the Sherwood number, containing ky, or ke, as a 
function of the Reynolds and Schmidt numbers, or their prod- 
uct, the Peclet number; alternatively, with a theoretic approach 
to be subjected to an empirical confirmation, such correlations 
can be derived from a physical—mathematical description. For 
instance, in Part 1, Appendix C, Eq. (C19) for the uniform 
velocity electrode is equivalent to the correlation 


Sh = 1.128 Re!/? Sc!/* = 1.128 Pe!/? (11) 
kL L 
wee pea A 
Da u 
L 
es peni a2) 
Dap Da 


2.3. Constant current analysis 


When the object of the analysis is more particularly focused 
on mass transfer resistances and their dependence on the flow 
conditions, the desired information can be more quickly and 
simply obtained by varying the inlet conditions. For instance, 
the measurements can be performed under the constant current 
restraint and by changing the velocity at the electrode through 
flow rate changes, so that the current density J becomes a con- 
stant, while the limit current Ji is a variable. 

For the identification of the diffusive resistances from volt- 
age measurements at constant current the linearisation of Eqs. 
(9) and (10) can be used [2], but it is often preferable to make a 
direct reference to the non-linear form of Eq. (7), allowing for a 
constant ką. However, it should be remembered that, for a com- 
plex electrode, ky, does not coincide with kę and, moreover, for 
measurements under different flow rates, the reference condition 
for the limit current is not evident or obvious. On the contrary, it 
is often preferable to choose another type if reference condition, 
i.e. a very high flow rate condition (that is a very low utilisa- 
tion factor), at which the transport resistances connected to flow 
dynamics become very low or even completely negligible. 

As usual, any refinement of the analysis passes through a 
more detailed description of the flow and composition fields 
along the electrode, that is the reliable calculation of the coeffi- 
cients A and Ay. Making reference to working conditions with 
constant current and changing utilisation factor, in Eq. (7) ka is 
a constant while k and ke = KLAL are variables, depending on 
ke. It is then preferable to write Eq. (7) as follows: 


1 1 1 1 1 (A—1)/A 
ko Aka ke ka ke ka 


(13) 


As already said, under these conditions it is not particularly 
useful to make reference to low currents, with kọ controlling 
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or high currents, with kL = kc/AL controlling, but it is better to 
refer to particularly favourable flow conditions (high flow rates, 
the current remaining unchanged), in which the transport coeffi- 
cient is high. Moreover, at high flow rates, the utilisation factor 
approaches zero and the averaging coefficient A approaches 
unity, while ke increases and 1/k, decreases as far as, in some 
cases, to become negligible. Under these reference conditions 
Eq. (13) attains its limit form 


—-— (14) 


and, by combining Eqs. (13) and (14), a useful identification tool 
is obtained 


koe _ 1 = Chalke) = (A= DIA) _ py) 
ko i= (ka/km) ea ore 
en 2 (15) 


by means of which the voltage loss is expressed, with respect to 
the reference, in terms of transport coefficients. 

Starting with polarisation data, Eq. (15) can be used in order 
to obtain the transport coefficient and its dependence on the inlet 
flow rate, on condition that an acceptable theoretical evaluation 
of the averaging coefficient A can be performed; inversely, Eq. 
(15) can be used to obtain information on the composition field 
starting with some kind of estimation of the transport coefficient. 
In the first section an example of the first type of approach will 
be discussed. 

Here a comment on small perturbations should be added: 
when the inlet flow rate is only slightly changed, Eq. (15) can 
be written in the simple differential form 


Pdi © — @ d In(ke) (16) 


on condition that ka < ke and, coherently, A ~ 1. In all the other 
cases it is probably better to make a direct reference to the finite 
formulation of Eq. (15). 

So, two difficulties can be underlined in using perturbative 
measurements; first of all they must be performed in terms 
of variations large enough to induce well measurable voltage 
effects in all instances; moreover, they are not always directly 
interpretable by neglecting the effects of the composition fields 
(A different from unity and depending on the flow rate). On 
the contrary, under the opposite condition, where the transport 
coefficient depends slightly on the flow rate (see the following 
section), the potential change is mainly due to differences in the 
flow field 


Bdn: ~ din(A) (17) 


and the evaluation of kę cannot be considered unrelated to an 
estimation of A. 


2.4. Analysis in terms of limit currents 


For each limit current condition, direct kų measurements lead 
univocally to the transport coefficient kę, on condition that the 


effect of the composition field, according to the second equa- 
tion in (6), is taken into account through the coefficient AL. 
Limit current data are particularly suitable because the estima- 
tion of A, is often simpler than that of the generic coefficient 
A, which is also a function of the ratio between the actual and 
the limit current. On the other hand, each characteristic curve 
gives one limit current value only, so that these data are relatively 
rare. 

However, it is worth underlining once more that the diffusive 
resistance cannot be automatically and directly determined from 
the limit current alone. 


3. An example of constant current analysis 


In Tables 1-3 and Fig. 1, some series of experimental data on 
an Ansaldo Fuel Cells molten carbonate fuel cell are reported. 

The experimental work was conducted at the Ansaldo lab- 
oratories; the experimental apparatus and test procedure have 
been described elsewhere [3]. Here two test series with a per- 
turbative nature are specifically discussed: for each series, in 
Tables 1 and 2, respectively, a central point is defined and the 
working range is reported in terms of different inlet flow rates, 
while the cell current is maintained rigorously constant at the 
central point value. The cell temperature and the humidification 
conditions are also maintained unchanged. Each series is articu- 
lated in two subseries and in each subseries the inlet flow rate to 
one electrode is varied while the flow rate to the other electrode 
is constant. Each experimental point refers to the stationary con- 


Table 1 
Ansaldo Fuel Cells: series I—test conditions 


Central point 
Anodic flow rates [N m? h7! x 107] (yao =0.225, A=Hz) 


Ho 7.43 
CO2 2.24 
N2 20.50 
H20 2.86 
Total 33.03 


Cathodic flow rates [N m? h7! x 10°] ao =0.161, A =02) 


Air 84.66 

CO2 14.31 

N2 6.23 

Total 105.20 
Cell dimensions (cm?) 10 x 10 
Cell current (A) 5.5 
Cell voltage (V) 0.814 
Anode 

H3 utilisation (%) 30.97 
Cathode 

CO; utilisation (%) 16.08 

Oz utilisation (%) 6.79 


Anodic tests 


Anodic flow rates 
Cathodic flow rates 


Cathodic tests 


Anodic flow rates 
Cathodic flow rates 


0.54an ep < qan < 1.5qancp 


aca = Ycacp = const. 


Jan = Jan cp = Const. 


0.54cacp <dea< 1.5qcacp 


350 P. Costa, B. Bosio / Journal of Power Sources 172 (2007) 346-357 


Table 2 
Ansaldo Fuel Cells: series H—test conditions 


Central point 
Anodic flow rates [N m? h7! x 107] (yao =0.727, A =H2) 


H2 7.80 

CO2 2.00 

No 0.00 

H20 0.93 

Total 10.73 
Cathodic flow rates [N m?h7! x 103] (Yao = 0.129, A =02) 

Air 91.50 

CO2 8.00 

No 46.50 

Total 146.00 
Cell dimensions (cm?) 10 x 10 
Cell current (A) 29 
Cell voltage (V) 0.768 
Anode 

Hp utilisation (%) 29.50 
Cathode 

CO, utilisation (%) 28.77 

Oz utilisation (%) 6.11 


Anodic tests 
Anodic flow rates 
Cathodic flow rates 


0.54an cp < Qan < 1.5qan cp 
Ica = 4cacp = const. 
Cathodic tests 


Anodic flow rates 
Cathodic flow rates 


qan = Yan cp = Const. 
0.54cacp <dea< 1.5qcacp 


Table 3 

Ansaldo Fuel Cells: series I and H—anodic results 

up, X° x’? V [mV] ye 

Series I, anodic—y4o = 0.225 
0.2010 0.250 0.209 828 0.735 
0.2382 0.297 0.249 825 0.708 
0.2817 0.352 0.297 818 0.648 
0.3097 0.388 0.329 814 0.617 
0.3440 0.433 0.368 809 0.579 
0.4425 0.564 0.485 785 0.429 
0.6186 0.812 0.716 723 0.198 

up, X x V [my] ye 

Series II, anodic—y4o = 0.727 
0.1967 0.328 0.204 789 0.698 
0.2269 0.376 0.237 783 0.648 
0.2682 0.442 0.282 773 0.571 
0.2950 0.484 0.312 768 0.536 
0.3278 0.535 0.349 759 0.479 
0.4215 0.680 0.460 729 0.330 
0.5901 0.938 0.677 645 0.116 


The bold data are related to the “central point” of the Series. 

àa X=2(A — 1)/A; A=—Yao — [A + yao Vu] ln(1 — u). 

8 X =2(A! — 1)/A'; A! = —(1/u) In — u). 

€ Y=exp[—A(n — n:)]. The assumption £= 1/2 is made and the reference 
potential is defined so that n =0 for A=1. 

d The choice of a unique external reference for both series I and II requires 
that the voltage data of series II are corrected by adding 34 mV, allowing for the 
Nernst concentration effect. 


40 


V-Vep 


©, @ anodic sub series 
-120 O, E cathodic sub series 


0.4 0.6 0.8 1 1.2 1.4 1.6 
d/d cp 


Fig. 1. Rough experimental results: voltage—flow rate trend for the four sub- 
series. x-Axis: ratios of the actual flow rate and the proper central point one; 
y-axis: differences between the actual cell voltage and the proper central point 
one. 


dition attained after each flow rate change. The results, in terms 
of cell voltage as a function of the inlet flow rate, are reported 
in Tables 3 and 4 and are shown as a whole in Fig. 1. 


3.1. Anodic tests 


The figure shows immediately that the cell is much more 
sensitive to the anodic flow rate changes; this is the reason why 
the analysis will be focused on the behaviour of the anode. 

An a priori evaluation of the order of magnitude of the diffu- 
sive resistances of the MCFCs built by Ansaldo Fuel Cells, and 
in particular the evaluation of the anodic resistances, has sug- 
gested that they are mainly localised in the condensed or still 
phases and, then, independent of flow dynamics (independent 
of the flow rate and the utilisation factor): resistances of this 
kind, independent of the flow pattern along the electrode, can 
be localised in the plate holes, the catalyst layer, and the liquid 
electrolyte. 


Table 4 
Ansaldo Fuel Cells: series I and II—cathodic results 


Series I, cathodic Series II, cathodic 


Uco, V [mV] uco, V [mv]* 
0.1072 808 0.1918 756 
0.1237 810 0.2213 756 
0.1462 810 0.2616 761 
0.1608 814 0.2877 768 
0.1787 809 0.3196 761 
0.2297 809 0.4109 760 
0.3214 801 0.5753 754 


The bold data are related to the “central point” of the Series. 

à The choice of a unique external reference for both series I and H requires 
that the voltage data of series II are corrected by adding 34 mV, allowing for the 
Nernst concentration effect. 
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By comparison, the resistances inside the homogeneous 
moving gas phase appear to be less important. So, in a first 
approximation, confirmed by Eq. (15), the assumption 


ke = kem = const. (18) 


can be introduced and the correlation 


A-1 
Y=1-CxX, Y=exp(—fy,), X =2 =| y 
ke 


CS. 
2(ke — ka) 


(19) 
is obtained. 

In such a way, the first equation in (19) foresees a linear 
relationship between the voltage variable Y and the composition 
variable X, the linearity constant C being directly connected to 
the transport coefficient. 

The coefficient A and, then, the variable X assumes a specific 
value at each experimental point, depending on the concentration 
distribution along the anodic side of the cell, especially for the 
hydrogen concentrations. These values are not easily accessible: 
a first attempt at a rough estimation should lead to assume the 
arithmetic mean Cam œ% Cao(1 — (u/2)), to which A © 1 + (u/2) 
and X7 u correspond; nevertheless, this first approximation is 
sufficiently correct only for very low utilisation factors. Fig. 2 
demonstrates that the identification of X with u does not allow a 
satisfying correlation of the experimental data. 

Now, if the single adduction channel of the cell is schematised 
as an ideal tubular reactor with an irreversible first-order reaction 
(see Part 1, Appendix C2 and Part 2, Appendix D) and, instead 
of the uniform flow rate assumption, the flow rate changes that 
are due to the hydrogen consumption (v=1) along the axial 
coordinate are explicitly taken into account (see Appendix E2), 
a reasonable estimation of the averaging coefficient A is 


1+ ya 
A =—YaAo | "| In(1 — u) (20) 
u 
0.8 
+ 
0.7 a? 
+ 
A te 
0.6 Š 
A 

n A 
= 0.5 4 
= + 
04 
a 
is N 

0.3 

0.2 + 

@ sub series | 
0 T T T 
0 0.2 0.4 0.6 0.8 


Fig. 2. The anodic polarisation as a function of the utilisation factor. 


0.8 


0.74 


0.64 


0.5 


0.4 


EXP(-Bn r) 


Fig. 3. The anodic polarisation as a function of the variable X¥=2[(A — 1)/A] 
with A =—yao — [1 + yao Vu] In — u). 


and only when both yao and u are very low can the variable X 
of Eq. (19) be confused with u: 
jig > 0, u —>0, ANI+5, Xxu (21) 
The data in Table 3 show that, with reference to the experi- 
mental conditions, the equations in (21) allow only a very rough 
estimation of the variable X. On the contrary, by using Eqs. (19) 
and (20), the estimations of A and X=2(A — 1)/A are substan- 
tially correct within the validity of the assumption made. It is 
useful to remember that these assumptions are: 


e isothermal and isobaric system; 

e anode effectively de-coupled from the cathode; 

e anodic flow pattern assimilable to that of an ideal tubular 
reactor; 

e irreversible kinetics, first order in respect to the hydrogen; 

e constant mass transfer coefficient. 


When X is calculated by means of Eqs. (19) and (20), the 
result of Fig. 3 is obtained, where the X — Y trend is shown to 
really correspond to the linear correlation of Eq. (19); moreover, 
a unique straight line allows the contemporary representation of 
the two subseries of data, which correspond to very different 
composition conditions. Only one warning is necessary, that is 
the choice of a unique external reference condition for both the 
subseries, so that the constant k, is really the same for all the data; 
choosing the inlet condition of the first subseries as reference it 
is necessary to add 34 mV to the measured voltage of the other 
subseries to allow for the less favourable inlet conditions. The 
linearity and the uniqueness of the obtained correlation are sig- 
nificant results. As said before, it is not possible to obtain similar 
results by calculating X through the oversimplified equations in 
(21) (see Fig. 2). 

In other words, using the equations in (19), with A from Eq. 
(20), ka is a known parameter, corresponding to the chosen con- 
stant current, A is a function of the utilisation factor, 7 is intended 
to be measured as a function of u. The transport coefficient ke 
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must be considered as an adjustable parameter, for which the 
most probable estimation is derived from the data correlation. 
In synthesis, this is a typical identification problem. When ke 
depends on flow dynamics, it can be expressed by equations 
like (11) in terms of Sherwood, Reynolds and Schmidt num- 
bers; then, the best-fit operation aims to determine the correlation 
exponents and especially that of the Reynolds number, in which 
the flow pattern effects are summarised. 

In the case under examination, from the slope of the linear 
correlation and the last equation in (19) the following simple 
result is obtained: 


ke © 2.lka © 5.2 x 1072 cms7! = const. (22) 


Thanks to the use of an unique external reference, Eq. (22) 
shows that the anodic diffusive resistances are significant, being 
about one half of the total, and really independent of the flow 
conditions. 

For this last purpose, it is useful to observe that the diffusive 
resistances can, in general, be split into two parts; this is equiv- 
alent to assuming that the transport mechanism consists of two 
paths in series, the first depending only on the geometry of the 
condensed or still phases and independent of the flow dynam- 
ics and the second regarding the moving gaseous phase and so 
depending on the flow dynamics: 

1 1 1 


eee +—, kp = const., keg = f(v) (23) 
ke kocp ke ° i 


A result such as that in Eq. (22) corresponds to the statement 
that the anodic diffusive resistance of an Ansaldo MCFC are 
prevalently of the first type, while the gas phase resistances are 
negligible: 


ke © kecp © keM, keg > kecp (24) 


In other words, the distributing plates close to the electrodes 
induce a gas flow pattern which is effective in minimising the 
diffusive path of the reagents from the moving bulk gas to the 
still electrode structure. 

It is remarkable that the structure of the coefficient kecp is 
not affected by the bulk flow leaving the anode (see Appendix 
E1, v=1). This circumstance, together with a rough evaluation 
of the resistances in the still gas phase (plate holes, gas filled 
part of the porous electrode) and in the liquid phase (wet part of 
the porous electrode), confirms that the last are controlling. 

In effect, the whole phenomenon is very articulated: the 
macroscopic flow fields are obviously relevant in the transport 
through the moving gas; the bulk flow effects should also be con- 
sidered as acting with the diffusive fluxes in the transport through 
the still gas; the crossing of the gas—liquid interface involves the 
solubility of the reagents; finally, the transport through the liquid 
phase should not be affected by the bulk flow, as this is certainly 
negligible due to dilution. A slightly more analytical discussion 
of the bulk flow effects is reported in Appendix E1. 

Allowing for the value assumed by the constant ką in the 
experimental conditions considered (kg ~% 2.5 x 107? cms~!), 
Eq. (22) leads to an estimation of the transport coefficient of 
the order of ke ~5 x 10-2 cms7!, which is at least an order of 


magnitude lower than that expected for a controlling gas phase: 
so from this point of view a controlling liquid phase appears to 
be very likely. 


3.2. Cathodic tests 


The analysis of the cathodic diffusive resistances is more 
complex and, at the same time, less interesting because the 
effects of cathodic flow rate changes are much more moder- 
ate. This fact can be related to one or more of the following 
considerations. 

In the first place, cathodic flow rates are decidedly greater 
than anodic ones, so that the flow field effects on mass trans- 
port are correspondingly less important. In the second place, 
there are some pertinent observations to be made on the com- 
parison of the diffusive losses, which are independent of the 
flow dynamics, with the intrinsic cathodic kinetics; it is known 
that this kinetics depends mainly on the oxygen concentration 
through a fractional reaction order (for instance, 0.5—0.7 [4]) and 
it is certainly slower (greater activation losses) than the anodic 
one; both these circumstances make the diffusive effects less 
appreciable. 

A third and probably conclusive effect regards the oxygen 
concentration changes during the reaction: the rather low utilisa- 
tion factor (about 6%) and a decreasing flow rate act to maintain 
the oxygen concentration almost constant (see Appendix E2). 

Therefore, while the anodic kinetics varies sensibly, the 
cathodic kinetics, at least under the experimental condition 
tested, should not vary by more than a few percentage points. 
The effects of the composition field on the longitudinal flow elec- 
trodes of a molten carbonate cell, the more relevant anodic ones 
and the much smaller cathodic ones, can be directly observed in 
Fig. El. 

The small voltage changes (less than 15 mV, that is about 
2%) observed in the cathodic subseries were in agreement with 
what was just said in order of magnitude. More particularly, the 
non-monotonic trends, provided they are significant, could be 
connected to little non-uniformity in the internal temperature 
field, which are not completely eliminated by the thermal reg- 
ulation, or, more probably, connected to coupling effects with 
the more sensitive anode. In other words, while the anodic per- 
turbations can be considered well de-coupled from the cathodic 
behaviour, the cathodic perturbations could modify the anodic 
behaviour more than they modify the cathodic behaviour. 

All that considered, the cathodic effects are too small and the 
acting phenomena too numerous and non-sufficiently controlled 
to risk a quantitative interpretation. On the contrary, the anodic 
effects are relevant and can be reliably attributed to the hydrogen 
concentration changes along the reaction coordinate. 


3.3. Results 


In short, from the above example of constant current analysis, 
the following conclusions can be drawn: 


e The mass transport phenomena at an MCFC anode can be 
determined by means of perturbative methods applied to the 
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inlet anodic flow rate, the cathodic conditions and the cell 
current delivered remaining unchanged. 

e The changes in the cell voltage can be directly related to the 
hydrogen concentration field along the anode; however, the 
effective mean concentration must be evaluated with some 
attention to the real flow conditions of the anodic gas, in 
agreement with chemical reactor theory fundamentals. 

e The anodic analysis results indicate that the anodic diffusive 
resistances are significant and can be considered comparable 
with the anodic activation resistances. 

e The mass transport phenomena and the composition field 
effects are much less evident at the cathode and, presumably, 
less important than the cathodic activation losses. 

e The diffusive losses appear to be negligible in comparison to 
the activation losses at the cathode and comparable to the acti- 
vation losses at the anode. As it is known that the activation 
losses at the anode are lower than those at the cathode, the 
concentration effects are, as a whole, responsible for a rela- 
tively small part of the total losses. However, as the anodic 
utilisation factor increases, these effects become more and 
more evident. 

e The diffusive losses within the moving gas are negligible, 
while the ones inside the liquid appear to be controlling. 
This, among other considerations, demonstrates the effective- 
ness of the distributing plates used in the Ansaldo Fuel Cells 
technology. 


4. An example in terms of limit currents 


Tables 5 and 6 report some experimental data, in terms of 
limit currents, which have been obtained under constant flow 
rate conditions and by varying the concentration of the anodic 
feed. These data have already been presented in a previous work 
[5], where a substantial proportionality between the limit current 
and the inlet hydrogen molar fraction (yao) was demonstrated. 
Those conclusions are now more precisely discussed. 

In agreement with what has been discussed in the previous 
section and allowing for the same assumptions on a longitudinal 


Table 5 
Ansaldo Fuel Cells: test condition 


Central point 
Anodic flow rates [mol s~! x 10°] 


H2 6.5° 
CO2 1.4 
N2 15.8 
H20 2.3 
Total 26.0 
Cathodic flow rates [N m? h7! x 10°] 

O2 15.1 
CO2 14.3 
N2 125.0 
Total 154.4 


Cell current Limit current 
Flow rates 


Hydrogen mol fraction 


Constant flow rates 
0.06 < yao < 0.70 


a H, and N3 flow rates are changed in a complementary way, so that the total 
anodic flow rate is kept constant. 


Table 6 

Ansaldo Fuel Cells: limit current results 

Yao Ji [mA cm~?] uL AL’ YAo/AL 
0.699 270 0.422 1.510 0.463 
0.583 260 0.487 1.588 0.367 
0.437 230 0.575 1.700 0.257 
0.350 205 0.640 1.804 0.194 
0.291 178 0.668 1.842 0.158 
0.250 153 0.669 1.812 0.138 
0.219 140 0.699 1.872 0.117 
0.194 120 0.676 1.796 0.108 
0.175 100 0.624 1.667 0.105 
0.146 100 0.748 1.972 0.074 
0.117 80 0.747 1.950 0.060 
0.097 68 0.766 1.980 0.049 
0.087 55 0.690 1.761 0.049 
0.079 55 0.761 1.950 0.041 
0.070 50 0.781 2.010 0.035 
0.062 40 0.705 1.777 0.035 


a Ar =[—yao —(1+yYao) Ind — ur X/uL]~'. 


flow anode, the relationship between the limit current (that is 
kL), the transport coefficient ke and the inlet composition yao 
should be expressed by 


k 

ie = —yYAoUL — (1+ yao) Ind — uL), 
k J 

u, =~ = — + (25) 
ke = nFyaoke 


Assuming that k and k¢ are constant, Eq. (25) can be written 
in the form 


n FkcYyAo const. yAo 
JL a = ’ 
AL AL 
1 1 — uL 
— = |—yao — (1 + Yao) In (26) 
AL uL 


in which the simple proportionality between JL and yAo is 
corrected in order to allow for the non-uniform concentration 
field. 


0 50 100 150 200 250 
JL 


Fig. 4. Effects of the inlet hydrogen mol fraction yAo and of the mean hydrogen 
mol fraction yAo/AL on the limit current. 
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Only in the case of very low utilisation factors can it be 
assumed that the corrective term Ay, connected to the concen- 
tration field is equal unity (uniform concentration) and can Jg 
be considered directly proportional to yao. In Fig. 4, the differ- 
ences between the two correlation methods are shown: on the 
one hand, both the JL — yao and the JL — yao/AL correlations 
indicate straight lines from the origin of the axes, but the two 
slopes are completely different. In other words, neglecting the 
factor A; results in a serious underestimation of the transport 
coefficient ke. 


5. Conclusions 


As concluded in Part 1 of this work, the analysis of the diffu- 
sive resistance in a finite-size electrode can lead to a non-linear 
problem at the laboratory scale, where temperature and volt- 
age can be assumed uniform on the cell plane. At least two 
interlinked factors have to be considered: 


e the local diffusive resistances set and 
e the concentration distribution on the electrode. 


The first can be determined only after estimating the second 
by means of an adequate physical-mathematical description. 
To do this it is particularly useful to refer to a close analogy 
between electrode kinetics and heterogeneous chemical kinetics 
with the reaction at the boundary and to use this analogy within 
the chemical reactor theory. 

This kind of analysis can be successfully applied to experi- 
mental fuel cell data. The consideration of data from Ansaldo 
MCEFCs makes it possible to say that 


e the diffusive resistance at the cathode has less relevant effects, 
while 

e the anodic resistances are comparable with the activation ones 
and 

e are independent of the inlet flow rate, that is 

e prevailingly localised in the liquid phase; 

e a superficial analysis of the concentration field effects, in 
terms of constant current and limit current, can lead to a gross 
overestimation of the mass transport coefficients. 
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Appendix D. Comparison of electrodes 


The reference to limit models, or ideal reactors, and their 
comparison are fundamental classical arguments of the chemical 
reactor theory. In terms of continuous systems, in particular, two 
ideal reactors are compared, the continuous stirred tank reactor 
(CSTR) and the plug flow reactor (PFR). The undoubted useful- 
ness of these methodological schemes in the chemical reactor 
analysis, both homogeneous and heterogeneous, suggests that 


they should be more frequently utilised in electrochemical reac- 
tor analysis too. In Part 1, Appendix C, the behaviour of the 
“mixed”, that is the CSTR-like, and the “longitudinal flow”, that 
is the PFR-like, electrodes have been discussed, especially with 
reference to their characteristics as complex electrodes; here fur- 
ther remarks will be made to clarify their differences from the 
electrochemical point of view. 

A comparison of the PFR and CSTR, when strictly considered 
as chemical reactors, highlights, as known, the superiority of the 
first, which presents a greater mean reaction rate. When dealing 
with electrodes, the comparison can be better appreciated if it 
is made in terms of voltage, the utilisation factors (that is the 
currents or the mean reaction velocities), the active surfaces per 
unit flow rate and the mass transfer resistances being the same. 
Then, for the two reactors, or electrodes, equal flow dynamic 
constant ks, transport constant ke (or at least their ratio) and the 
utilisation factor u are assumed. 

A first difference in favour of the longitudinal flow reactor 
regards the limit utilisation factors, which are, for the two cases, 
respectively 


ke 
PFR u= l- exp| -> 
ke 
_ 1 
1 + (ke/ke) 
the first being greater than the second: the comparison at equal u 
cannot be extended to the range of higher values only accessible 


to the PFR. 
In the range where the comparison is possible, it is 


(D1) 
CSTR uL 


1 ( =) : 5 ee (D2) 
u = ex = x = 
PUE IFD è Kk k 
or 
k 1 k 
ii = ~ hA) 
o nud —u c 
kf l—u kẹ (D3) 
CSTR “£ = esis 


ky u ke 
At equal k¢/k, and u, the longitudinal flow, thanks to a bet- 
ter concentration distribution, requires a lower intrinsic kinetic 
constant ko, which means lower activation losses. The logarithm 
of the ratio of the two ky values can be directly connected to the 
differences between the activation losses of the two electrodes 


BAn=I1n |- 


= 4 1 
—In —-—|, u < ———_ 
1+ (ke/ke) 


(D4) 


The trend of this difference, which tends to diverge when 
the mixed electrode is approaching its limit current condition, 
is reported in Fig. D1 for different values of the coefficient ratio 
kelke. 

From another anymore relevant point of view, at least here, 
the comparison of electrodes presenting different flow charac- 
teristics involves the interpretation of the experimental data and, 
in particular, affects the analysis of the diffusive resistances: the 
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1 k,/k, = 20 


0.4 0.5 0.6 0.7 0.8 0.9 1 
u 


Fig. D1. Comparison of a longitudinal flow electrode (PFR) and a mixed elec- 
trode (CSTR) in terms of voltage losses as a function of the utilisation factor 
and for different values of ke/k¢. The difference An between the CSTR losses 
and the PFR ones is reported: each unity on the ordinate scale corresponds to 
about 40 mV (ke/kf = 20, utcstR = 0.95; ke/k¢ = 10, uLcstR =0.91). 


characterisation of one electrode in terms of exchange current 
and limit current does not allow a univocal estimation of the 
transport coefficient, unless the electrode is a simple one. In the 
simple cases it is straightforwardly 


ae : ke =k (D5) 
ko = ka kE’ c = KL 
and the transport coefficient is directly connected to the limit 
current. 
A continuous mixed electrode (CSTR) with given exchange 
and limit currents behaves just like simple electrode (see Fig. D2) 
1 1 1 1 1 1 
= , = (D6) 
ko ka k ke ke ke 
but its limit current also depends on the flow constant kf and, 
inversely, the limit current is not directly interpretable in terms 
of a transport coefficient. 


0 0.2 04 06 08 1 1.2 
We 


Fig. D2. A comparison of the electrodes in terms of voltage at different cur- 
rents, the exchange and the limit current being unchanged. The electrode voltage 
reported is in mV. 


4 


simple 


0.1 0.2 0.3 0.4 0.5 0.6 0.7 
k/ks 


Fig. D3. Comparison of electrodes in the identification problem: the estimation 
of the transport coefficient at equal exchange and limit currents. 


The longitudinal flow electrode (PFR), at equal exchange and 
limit current, also shows qualitatively similar performances (see 
Fig. D2) 

1 1 i 1 
ko ke In. — (ka/k) ke In. — (kL/ kò)’ 


ke (D7) 


ll 
| 
> 
SS 
— 
5 
—, 
_ 
| 
ye 
Sau 


and the limit current is again affected by the flow constant, but 
in a different manner from the mixed electrode and with a low 
current behaviour which risks being attributed to a different 
exchange current. 

As an example, in Fig. D3 the ratio kc/kų between the trans- 
port coefficient and limit coefficient of the three electrodes 
considered (simple, CSTR, PFR) are reported as functions of 
the ratio k/ke between the limit and flow coefficients. It is quite 
evident that the identification of a transport coefficient may be 
completely misinterpreted if these flow effects are not correctly 
taken into account. 


Appendix E. Bulk flow effects 


In fuel cell electrodes the electrochemical reactions are usu- 
ally unbalanced in terms of mol numbers and in terms of mass. 
This fact might add significant non-linearity effects to the local 
kinetics as well as to the concentration distribution on the cell. 


El. Local effects 


As is well known, the molar fluxes of the reagents from the 
bulk gas to the electrode surface and those of the products in 
the opposite direction contain two terms. The first, of diffusive 
nature, can be expressed through the Fick equation, while the 
second is a bulk flow term, which is due to the mean molar veloc- 
ity of the fluid, that is the total flux, the sum of all component 
fluxes: 


Ni = Nyy; — CrD : Nr=) Ni (El) 
i = iNT yi TVA dz ; T= i 
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On the other hand, the reaction stoichiometry states that 
Nj=vjNa/va, so, for the key reagent migrating towards the 
electrode, Eq. (E1) can be written 


dya/dz ) e ovi E2) 


1 + vya VA 


Na = —CryDa ( 


and, by integrating on a thickness s between ya; and yap, the 
flux expression is obtained 


1+ 2a] (= 
1 + vya2 


7 (=>) i we (E3) 
S 


which is non-linear in the composition difference. When v 4 0, 
the bulk flow effect is negligible only for very dilute solutions. 
Table E1 reports the values for v for MCFC, PEMFC and SOFC 
electrodes. For v>0 the bulk flow contrasts with the motion 
towards the electrode, so that the key reagent has a flux lower 
than the one associated with a simple diffusive path (v =0); for 
v <0 the effects are the opposite. 

For molten carbonate fuel cells and, in particular, under the 
experimental conditions reported in Tables 1 and 2, these local 
effects, especially the anodic ones, should be present and appre- 
ciable. If, on the contrary, similar effects have not been detected, 
a likely interpretation could be found in the role of the liq- 
uid phase. If two phases, gas and liquid, are considered, their 
separation interface can be assumed to be close to equilibrium 


Na = ln | ) = (yal — ya2) 


vs 


YAi = MAXAi (E4) 


Moreover, the molar flux of the key reagent is the same in each 
phase 


CgD 
Na = (YA — Yai) ( Pee) fe, 


and equal to the reaction rate per unit surface 
Na = kCLPXAs (E7) 


Under such conditions, the transport phenomena are con- 
trolled by one or the other phase, according to the value assumed 
by the ratio 


ma foCcDac/sc 


(E8) 
fipCLpDacp/SLP 


As the gaseous species A has low solubility (ma > 1) and 
CoDac > CLP Darp, comparable resistances in the two phases 
require a liquid thickness much lower than the gas one (for 
instance, SL < 10-45), so that it is quite probable that the liquid 
is the controlling phase. In this case, as ya; = maxai © ya and 
putting yas =MAXas, the flux is 


CipDap 
Na © (YA — Yas) ( ) fir (E9) 
MASLP 
where the global transport coefficient 
Cip D 
ka LPYALP (E10) 
MASLP 


depends on the transport properties and geometry of the control- 
ling phase and on the solubility coefficient ma, while the bulk 
flow effects become quite negligible because of the high dilution 
of A in the liquid phase 


Xai, XAs XK 1, fip © 1 (E11) 


Therefore, the absence of evident bulk flow effects in local 
diffusive phenomena is not surprising; on the contrary, this fact 
can be considered as a strong indication in favour of the liquid 
controlling assumption. 


E2. Effects on the composition field 


-h v(ya + yai) E5 i 
f= |l- 2 (E5) The molar fluxes unbalance shows a more appreciable effect 
on the composition of the bulk gas flowing on the electrode. By 
N Cup DALP expressing the molar flow rates of the various components as 
A = (Ai — XAs) SL Sip, well as the total molar flow rate in terms of the utilisation factor, 
: . Fsotvi 
fir = h 7 rataa (6) F=Fo- 2,  Fr= Fro + Faouv (E12) 
2 A 
Table E1 
Molar and mass unbalance for electrode reactions 
Electrode In Out A v= -X vi/va Vm = -X v;MilvaMa 
MCFC 
Anode H2 CO2, H20 H2 1 30 
Cathode 2C02, O2 - O2 —3 —3.75 
Cathode 2C02, O2 — CO2 —1.5 —1.87 
PEMFC 
Anode Ho - H2 —1 —1 
Cathode O2 2H20 O2 1 0.125 
SOFC 
Anode H2 H20 H2 
Cathode O2 — O2 —1 —1 
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the molar fraction of the key reagent is 


_ l-u 
YA = YAo 1+ yAck 


and the dependence on the molar unbalance is evident. 
In general, the dependence of kinetics on the key concentra- 
tion can be expressed in terms of an apparent order y 


(E13) 


Na =r = koCry%, = koCry%, O<a<1, a<y<l1 


(E14) 


which, in comparison with the true order œ, summarises the 
importance of the transport phenomena (y =a for kinetics con- 
trolling; y=1 for transport controlling). By combining Eqs. 
(E13) and (E14) the dependence of the kinetics on the bulk gas 
composition (ya,) becomes 


1—u)]|” 
cee Feared A E i 70 tad 
1+ vyaou 
(E15) 
1 dr 
? C) (=) TEE (E16) 
r du u=0 


For an MCFC the effect is much more important at the anode 
(v=1 and, with reference to Tables 1 and 2, yao = 0.2-0.8, 
u<0.3), where the corrective term on the right hand side of 
Eq. (E15) leads to reaction rate variations up to 60% in com- 
parison to the initial values; the effect is much slighter at the 
cathode (v=—3, yao ¥ 0.16, u% 0.06), where the decreasing 
flow rate contributes to maintaining the reaction rate substan- 
tially unchanged, with differences of less than about 3% of the 
initial value. 

The corrections to the local kinetics (Eq. (E15)) due to the 
unbalancing parameter v 40 also obviously affect the averag- 
ing operations and the values obtained for the parameter A, as 
defined by Eq. (5). For instance, in those cases where linear 
kinetics (w=y= 1) and a longitudinal flow can be assumed, 
the differential balance equation along the adimensional axial 
coordinate ¢ 

du’ Kd — u’) 


k— = l 
de 1+ vyao 


c=0, uv =0, c= 1, w =u 


(E17) 

leads to the integral 
K 
gy, 7 AO (1 + vyao) Ind — u), ka = ku, 

f 

A= > (E18) 

E 
and, finally to 
ln(1 — u) 

A = —vydo — (1 + vyao) ——— (E19) 


3 


CAo/CAm 


A= 


Fig. E1. Effects if the inlet composition and the flow rate changes on the MCFC 
longitudinal flow (PFR) electrodes (anode: v= 1; cathode: v = —3). 


This equation, with v= 1 is the one used in Sections 3 and 
4 for an MCFC anode. For a very diluted system as well as for 
balanced reactions (v =0) the same Eq. (E19) becomes 


Yao > 0, AX — (G) In(1 — u) (E20) 


u 
which is the one used in Part 1, Appendix C for describing the 
longitudinal constant flow electrode. Moreover, for very low 
utilisation factors, Eq. (E19) can be rewritten as 


u> 0, A~ 1+ (1+ vyao) (E21) 
while the equality of A with the mean utilisation 
vyao > 0, u—> 0, ANI+5 (E22) 


is correct only when both v, yao and u are very small. 

The effects of the molar unbalance on the averaging coeffi- 
cient A of the two longitudinal flow electrodes of an MCFC are 
reported in Fig. E1. The greater relevance of the anodic effect 
is evident. The cathodic effect is much smaller and may even 
change its sign at an increase in the oxygen content (yao > 0.33). 

In practice, the unbalance effects in MCFC are even lower 
than those of Fig. El because an intrinsic kinetics where O2 
is the key reagent (v = —3) interacts with transport phenomena 
where CO3 is likely to be controlling (v = — 1.5). 
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